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\^ ■ Abstract 



■ In this lecture at a school for condensed matter physicists, I be- 

I gin with basic concepts and tools for investigating phase transitions 

in quantum field theory. The very different roles of global and gauge 
. symmetries in phase transitions will be elucidated. Among the im- 

0^ I portant applications of the basic theory the thermodynamics of the 

electroweak transition is treated in detail, and the implications for 
O ■ baryogenesis will be discussed. The status of lattice simulations of 

the transition from a high temperature quark-gluon plasma to con- 
Qh' fined hadronic matter with spontaneously broken chiral symmetry is 

Q ■ also briefly reviewed. The various changes of the free energy density 

in cosmological phase transitions reinforce the cosmological constant 
^ I problem. In the last part of the lecture I shall address this profound 

mystery of present day physics. 

X. 

1 Introduction 

During the evolution of the very early Universe there have been at least 
two phase transitions. The electroweak theory predicts that at about 100 
GeV there was a transition from a 'symmetric' high temperature phase with 
'massless' gauge bosons to the Higgs phase, in which the SU (2) x f/(l) gauge 
symmetry is 'spontaneously broken' and all the masses in the model are 
generated. A detailed understanding of the nature and dynamics of this 
transition is a very difficult task, and a lot of quantitative analytical and 

*Invited lecture at the 3'"'^ Summer School on Condensed Matter Research, 7-14 August 
2004, Zuoz, Switzerland. 
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numerical studies have been performed over the years. One of the motivations 
of these studies has been that non-perturbative effects might have changed 
the baryon asymmetry in the phase transition. 

Quantum chromodynamics (QCD) predicts that at about 200 MeV there 
was a transition from a quark-gluon plasma to a confinement phase with no 
free quarks and gluons. At about the same energy we expect that the global 
chiral symmetry of QCD with massless fermions is spontaneously broken 
by the formation of a quark pair condensate. These transitions have been 
studied extensively in the framework of lattice gauge theory. It is unclear 
whether they have interesting cosmological implications. 

In the framework of unified extensions of the standard model other poten- 
tially very interesting phase transitions are expected. In some cases these are 
accompanied by the formation of topological defects (domain walls, strings, 
monopoles and textures). This can be cosmologically interesting or danger- 
ous. For example, in a second order phase transition of a grand unified theory 
(GUT) one expects generically that a completely unacceptable abundance of 
magnetic monopoles would be formed. This so-called monopole problem was 
one of the early motivations of inflationary cosmology. The first inflationary 
model by A.Guth ('old inflation') was based on a first order phase transition, 
and made use of the enormous energy density of the 'false vacuum' for the 
unbroken phase. 

These introductory remarks indicate the potential importance of phase 
transitions^ in cosmology. In this evening lecture I can, of course, treat only 
some selected aspects of the subject. Since this is a school for condensed 
matter physicists, 1 shall begin with foundational material for investigating 
phase transitions in quantum field theory. A central tool of most studies is 
the path integral representation of the partition function. In this connection 
I want to make some critical remarks about the widespread use of non-convex 
effective potentials, especially in the older literature, that disturbs me very 
much. 1 also want to explain why the roles of global and gauge symmetries in 
phase transitions are very different. In particular, a gauge symmetry can not 
be spontaneously broken (Elitzur theorem). As an important application of 
the basic theory, I shall review existing work on the thermodynamics of the 
electroweak phase transition, and its role in baryogenesis. We shall continue 
with a few remarks about the current status of lattice simulations of the QCD 
phase transitions. The subject of cosmic topological defects that might have 
been formed in phase transitions will only be touched. 

shall sometimes use the term phase transition in a loose sense, even if a detailed 
analysis may later show that the transition is a crossover, rather than a bona fide phase 
transition. 
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The various changes of the free energy density in cosmological phase tran- 
sitions reinforces the cosmological constant problem, a deep mystery of present 
day physics. In the last part of my talk I shall address this topical issue, and 
then summarize the current evidence for a dominant component of what peo- 
ple (unfortunately) call Dark Energy in our Universe. Whether this is a tiny 
remainder of an originally huge vacuum energy in the very early Universe 
remains open, but observations may eventually settle this issue. 

2 Phase Transitions in Quantum 
Field Theory 

Let us first consider a quantum field theory without local symmetries. The 
partition function can (formally) be represented as an Euclidean path (func- 
tional) integral. If S[(j)] denotes the Euclidean action functional of a (multi- 
component) scalar field 0, say, and if we include the coupling to an external 
source J(x), this reads^ 

Z[J]= f exp{-5[</.] + (</., J)} [#], (1) 

where the integration extends over fields which are periodic in Euclidean time 
r, with periodicity {3 = 1/kT. For textbook treatments see P and [2]. 

Several remarks on this formula are in order. In one dimension, i.e. for 
quantum mechanics, it has - suitably interpreted - a precise meaning and can 
be rigorously derived for a large class of interactions. In four dimensions this 
is rarely the case without an ultraviolet cutoff. An exception are free fields 
(Gaussian probability measures on spaces of distributions). In perturbation 
theory, renormalization theory leads - for renormalizable field theories - or- 
der by order to meaningful expressions. In some cases, for example for QCD, 
we expect that the continuum limit of lattice regularizations exists and is 
physically sound. This is, however, not the case even for simple scalar mod- 
els (triviality theorem |3]). It is remarkable that this Euclidean formulation 
in terms of path integrals establishes such a close connection with classical (!) 
statistical mechanics (models of magnetism). For the study of the thermo- 
dynamics of quantum field theories and other issues one can work entirely in 
the Euclidean framework. An explicit continuation back to Minkowski space 
is not needed. 

^We consider here only the canonical ensemble. For grand canonical ensembles one has 
to replace in the action 84 by — /x, where /Lt is the chemical potential of the field. 
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An important quantity is (— /3) times the free Helmholz energy (Schwinger 
functional) : 

W[J\^\^Z[Ji (2) 

that generates the truncated correlation functions. The effective action 
(Gibbs potential) F is the generalized Legendre (Legendre-Fenchel) trans- 
form of W^[J]: 

r[$] = sup{($, J) - W[J]}. (3) 
J 

Expectation values of observables A[(f)] are given by 

{A) J = Z[J]-' J exp{-5[0] + (0, J)} [dcp]. (4) 

W[J] is convex in J. For a lattice regularization this follows simply from 
Holder's inequahty (see footnote below). Formally, this is imphed by 



j|^l^/Wd^.A = ((W)-(W)))')>o, 



for every real test function. As a Lcgcndrc transform, r[$] is also convex. 
For a homogeneous source J the effective potential Ve//($) is, by definition, 
the free Gibbs energy density. 

Consider, for illustration, the simple example of an 'action' without a 
kinetic energy term 

S[(l)] = = J U{(t)) d^x. 

For T = and a constant source J, the free energy density Wq = 1^0/^4 
in a finite region of volume of a lattice regularization is given by a one- 
dimensional integral: 

Wo{J) = 7^ In / e^'^'^'^-^^^'^^dif. (5) 

In the limit Q4 ^ 00 the Laplace approximation becomes exact, hence 

Wo^L{U), (6) 

where L denotes the Legendre-Fenchel transformation. This agrees (see Ap- 
pendix A) with L of the convex hull of U, and since L{L{U))= convex hull of 
U, we obtain the result that the Gibbs free energy density is for a constant 
$ equal to the convex hull of U. In Fig. 1 we show for a non-convex U 
the transforms L{U) and L{L{U)). L{U) is still the free energy density for 
T > 0. (Note that for T > the volume is /3 times the three- volume VL^.) 
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Figure 1: Graphs of a non-convex potential U, together with its first and 
second Legendre-Fenchel transforms. 

Phase transitions 

In Figs. 2 and 3 we show the typical shapes of effective potentials for first 
and second order phase transitions. Their Legendre-Fenchel transformations 
are the Helmholtz free energy densities. The reader is invited to sketch these. 
For instance, the slope of w{J) corresponding to Fig. 2 develops for T | Tc 
and J = a jump (discontinuity in the 'magnetization'), while for Fig. 3 the 
'magnetization' vanishes at T > Tc and starts to grow below T^. 

Footnote: Convexity of W[J] 

The convexity of W[J] holds even in a finite volume, and not only in the 
thermodynamic limit. (Actually, W is then even strictly convex). The lattice 
regularized partition sum is then an integral of the form 



where cx.={ai, ...,ak) and a • f{x) = J^o^ifii^)- We claim that for any A 
with < A < 1 




\nZ(^Xcx + (1 - A)/3) < X\nZ{a) + (1 - A) \nZ{f3). 



This follows from the following application of Holder's inequality 



FG\U< IIFIIpIIGII,, - + - = 1 (p>1). 



p q 



for F = e^f, G = e^^-^^^, p = 1/A, q = 1/(1 - A), e L\^). 
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Vpff{<3>,T) i 




Figure 2: Effective potential for a first order phase transition. At the critical 
temperature the minimum is degenerate (horizontal line corresponds to a 
phase mixture). The order parameter makes a jump between the pure phases. 




Figure 3: Effective potential for a second order phase transition. For the 
critical temperature the second derivative vanishes at the minimum. The 
horizontal straight line below shrinks to a point at T^. 
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Usually, the convexity of In Z is established by considering the quadratic 
form 

g^6«, - (A»'^> > 0, 

where := YLi fi^i^ ^-^d the angular bracket denotes the expectation value 
with respect to the probability measure Z(q:)~^ exp^^o; • f^dji. 

This argument shows that In Z is convex as a function of any parameter 
that appears linearly in the exponential under the integral sign in Z. 

2.1 Mean field approximation 

Countless expressions and plots in books and articles violate the basic con- 
vexity property, usually as a result of some approximate calculations that are 
analytically extended to regions where they are no more valid. Even in some 
really good books that happens already for the mean field approximation. 
For illustration let me explain what is not done properly. 

Consider a one component neutral scalar field with Euclidean action in 
D dimensions 

= j [^(90)' + d'^x ^ r[0] + (7) 

To have well-defined expressions interpret all the formulae in the lattice regu- 
larization for finite (hypercube) subsets A of the lattice. (Then T[0] contains 
nearest neighbor interactions.) Beside the action S we consider an 'approxi- 
mate' action 5*0 and the corresponding probability measures 

d^xM = Z;'e-^'[d(t>l Zo^J e-^°[d<f>] (8) 

for all finite subsets A. The corresponding expectation values are denoted by 
{A)q. The starting point is Jensen's inequality 

(expA)o>exp((A)o) (9) 

for A — —S + {(f), J) + Sq, and for 5*0 we choose 

So ^ -{<!>, K)+U[<f>], (10) 

where K{x) is another external source. This gives, if Wq = In Zq, 

W[J] >Wo + {So-S+ (0, J))o = Wo[K] - {T[<P])o + ((J - K), (0)o). (11) 

For (f) :— (0)o we have 
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Since {(j){x)(f){y))o = {(p)o{(p{y))o, we find 



W[J] > W[J, K] := Wo[K] - r[0)] + ((J -K),4>). (13) 

The mean field free energy Wmf[J] is the supremum in K of the Landau 
functional W[J,K], often also called the Landau free energy. Since lVo[-ft'] 
is convex, K and are according to (12) in one to one correspondence. Hence 
we can write 

W^mf[J] = sup{(0, J) - (FoM + T[0])}, (14) 
</> 

where Fq is the Legendre transform of Wq. This equation shows that W^m/[<^] 
is the Legendre- Fenchel transform of Fo[0] + T[0], and is therefore convex. 
Furthermore (see Appendix A), the mean free Gibbs potential is given by 

Tmfi^] = convex hull {Fo[$] + T[$]}. (15) 

It is crucial that this is convex. Most authors work with the standard Leg- 
endre transform and end up for F^/ with the functional in the curly bracket, 
which is in general not convex. (For a typical example, see Figs. 24.1, 2 in 
jl].) Note also that the Landau functional, expressed as a function of J and 
(f), is given by the curly bracket in (14) 

>V[J,0] = (J,0)-(Fo[0]+r[0]). (16) 

For J = its convex hull is according to (15) equal to — F^/- We emphasize 
that the Landau functional is in general not convex. 

As an exercise, one may work out everything for the Ising model. 



2.2 One- loop approximation 

The loop expansion of the Gibbs free energy F[$] is a much used systematic 
approximation scheme. We recall here only the zero- and one-loop terms. 
These are obtained with the steepest decent method. 

The saddle points 0o of the exponent in (1) satisfy the classical equation 
of motion 

Let us first assume that there is only one such saddle point for a given J. 
This is, for instance, the case for an action of the form (7) for a convex 
potential t/(0), because the previous equation then becomes 

-A(/)o + f/'(0o) = J, 
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which has at most one solution if J{x) vanishes at infinity and (pQ is assumed 
to go to zero at infinity. In the expansion 

- (0, J) = S(<f>o) - (00, J) + l(^l^,A^|^) + 
with ip := 4> — (pQ, and 

Aix,y) = ^ s , (17) 

S(l){x)S(l){y) 

we leave out the higher orders and arrive at 

Z[J] ~ exp{(0, J) - S[<Po]} J e--^^^'^^\diP]. (18) 

We interpret all the formulas in the lattice regularization (for a finite sublat- 
tice). The leading order of VFfJ] is given by the Laplace approximation 

Wo[J] = (00, J) - S[<Po] = {LS)[J]. (19) 
Therefore, the leading order To of F is given by 

Tom = sm (20) 

('Landau approximation'). The one- loop correction is determined by the 
Gaussian integral 

|e-^(^'^'^)[dV'] = (det^)"^ 

For W[J] we arrive at 

W[J]^Wo[J]+Wi[J] + ... , W^i[J] = -^lndet(^). (21) 

For a given J the conjugate variable $ differs from the classical solution 0o 
by a first order quantity. Up to second order we can thus replace 0o by $ in 
the expression (21) for Wi. Therefore, the Legendre transformation (3) 

F[$] = sup{($, J) - LS[J] - WiiMJ)] + •••} 
J 

is given by 

F[$]=Fo[$] + Fi[$] + ... , Fi[$] = ^lndet(^), (22) 
where A is now the expression (17) at = 
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If we set in (7) ^7(0) = |m + -P(0), where P{(j)) is a polynomial (of 
degree < 4 for renormalizable interactions), then the 1-loop contribution of 
r($) is given by the following expression (up to $ and T independent terms) 



r.(ft*) = ilndet{^ 



(23) 



The temperature dependence is contained in the Laplacian (periodicity in 
the Euclidean time direction) . The functional determinant can be worked out 
formally in the continuum theory or from the lattice regularization, for which 
it is well-defined and can easily be computed using Fourier transformation 
on a finite lattice. The details are given in Appendix B, where we also show 
that the continuum limit of the temperature dependent part Vi (/?, $) of the 
effective potential (free energy density) exists, and is given by 



where 



fl{k) = Vk^ + M2, :=m^ + P"($). 



(24) 



(25) 



The reader is invited to generalize this to the grand canonical ensemble of a 
complex scalar field. 

The T = contribution requires renormalizations. As a starting point 
we use the lattice regularization, derived in Appendix B. According to (63) 
we obtain in the thermodynamic limit for an x-independent $ 



lim — In det 
i Z4 



In 



M2 + -^(l-cosA;,a)], (26) 



where the integral now extends over the four-dimensional Brillouin zone 
BZ = [— ^, ^]^. In what follows it is simpler to use a spherical cut-off 
^2 ^ ^^2. In the continuum limit a — > we find for the regularized one- 
loop effective potential Vi{^) at T = 0, up to independent contributions. 



2 JA;2<A2 



In 



k^ + V^'{<^) 



d^k 



(27) 



(Vo = [/), where we have introduced an arbitrary mass scale ji (leading to a 
change of Vi which is independent of $) . This can be written as 
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with 



\n{x + y)x dx. 



We note that the third derivative of / is given by 



fiy) 



{x + yf 



dx 



and converges for A 
polynomial in y. So, 



oo to 1/y. Hence, f{y) — |y^ln(y) + quadratic 



647r2 



+ quadratic polynomial in Vq / /i^ 



The second order polynomial can be compensated by a renormalization of 
the couphng constants in Vq — U. Note that we also have to renormahze a 
^-independent vacuum energy. This will become a crucial issue when gravity 
is included (see Sect. 4). 

After renormalization the effective potential for T = is, up to one loop, 
given by 

1 _ _ r, _ 

(28) 



1 



647r2 

Remark. A change of the scale /x^, /i^ — > can be absorbed by a 

change of the coupling constants in Vo($;p(t)), where the parameters p{t) 
satisfy, up to higher orders, the differential equation 



Vi"(*:p(i)) 



l2 



If 



TTt^ f Q 

we obtain for i = the renormalization group equation 

dg 3 



dt 327r2 



g , etc. 



(29) 



Convex potentials 



For a convex potential like [/(0) = \m'^(j)'^ + ^0"', M^(<l>) = + 3A<I>^ is 
positive, and the one-loop effective potential is convex. The minimum occurs 
at $ = 0, and in this symmetric phase the T-dependent effective mass is 



m [T) —m + 



(30) 
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The T-dependent term can be read off from (24) 



= 3A 



$=0 



d^k 1 1 
(27r)3^e/^^ -1' 



UJ 



= VW+rr^. (31) 



For high temperatures we have 



4 L 



(32) 



For m = the scalar field acquires a mass \/AT/2 which plays the role of an 
infrared cutoff. 



Problems with non-convex potentials 

For a non-convex potential U (with the 'wrong' sign of the mass term) we 
can not simply use the previous formulae, since M^($) becomes negative for 
small $, whence Vi{j3, $) in (24) becomes complex. This means that the one- 
loop potential is not valid for small arguments, that are of particular interest 
in studies of the early universe. What should one do in this situation? We 
shall address this question below. 

For a rough estimate of the critical temperature for the transition to the 
symmetric phase we may naively ignore this and use the leading corrections 
in (32): 

m\n) = ~ + ^XT^. (33) 
This gives (Kizhnit, Linde, Weinberg) 

T2^t^ (m2<0). (34) 
A 

For weak coupling this implies that ^ — m^, and we can therefore expect 
that (34) is a good estimate. 

The unrenormalized parameter in this equation can be replaced by the 
mass of the physical particle (the "cr meson"). In lowest order, = — 2m^, 
hence 

T^--^. (35) 

We can, however, not ignore the encountered problems with the one-loop 
effective potential. It is quite clear that the saddle point approximation 
can only hold near the minima. If there are, for example, two similarly deep 
minima, as is expected close to the critical temperature, one has to take both 
into account. As a recipe one might take the convex hull of the sum of the 
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saddle point approximations for the two minima. But this is just a reasonable 
recipe. Another approach was studied in [S], where the standard definition 
of the effective potential is modified by coupling the source to a quadratic 
polynomial of the scalar field. The new potential is not necessarily convex 
as a function of the average value of 0, but is convex as a function of the 
composite field since this is the conjugate variable. The modified potential 
closely tracks the usual one where the latter is real. Moreover, at finite 
temperatures it displays a local minimum at $ = 0. It is also worthwhile to 
note that the resulting critical temperature in the one-loop approximation 
confirms the rough estimate (34). 

2.3 The constraint potential 

We consider now the free energy for a given spatial average of the order 
parameter. For simplicity, we do this for a scalar field theory, since general- 
izations are obvious. As before, the formulae below are well-defined in the 
lattice regularization. For a finite region A we sample the field configurations 
according to their mean field 

The probability distribution of for a homogeneous source J is given by 

= ^^exp(|A|[J0-K(0)])#. 

The constraint potential Vc (which depends on A) is given in terms of 
the integral in the last equation as 

^-WVM = I e-^AM5(0 - ^ / d^x) [c/0]a. (36) 

Since dP is a probability measure, we have 

Za(J) = J el^lt-'^-^^^^^lfi^. (37) 

In other words, W^a(j) is obtained from by a Laplace transformation. For 
finite regions Vc is usually not convex. 
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In the infinite volume limit the Laplace approximation again becomes 
exact if the saddle point of — J0 is unique in this limit. Hence 



which implies that the effective potential is equal to the convex hull of the 
thermodynamic limit of the constraint potential. Since it can be shown 
pn] that the latter is convex (certainly if the saddle point is unique, but 
presumably in general), the two potentials must be the same in this case. 

Note that for J = the probability density to find the system in a 
state of 'magnetization' (p is given by 



The constraint effective potential has thus a very intuitive meaning. It has 
often been studied via lattice simulations. For an example of an interesting 
recent application (vacuum stability), see [TT] . 

2.4 Effective potential for gauge theories 

For gauge theories it is advisable to construct a gauge invariant effective 
potential, that has a physical interpretation in terms of a gauge invariant 
order parameter. In this context it has to be stressed that for lattice regu- 
larizations of locally symmetric field theories any local quantity which is not 
gauge invariant, and has no gauge invariant component (for precisions, see 
Appendix C), such as a Higgs field, has a vanishing expectation value at any 
temperature. In this sense gauge invariance can not be spontaneously broken 
(Elitzur theorem [^). In perturbation theory, where one expands about a 
Gaussian measure, this is not so because gauge invariance is explicitly broken 
by gauge fixing. In the lattice formulation one does, however, not have to 
fix a gauge. (In practice, it may be useful to fix the gauge if the ensemble of 
configurations is not very large.) 

Lattice formulation of gauge theories 

The lattice regularization of gauge theories (Wilson) is of crucial importance 
in the strong (confined) field regime, because there is no other way to perform 
numerical calculations. But also from the conceptual point of view, lattice 
gauge theories are very interesting. For systematic presentations, see [Zj and 




1 





0. 



14 



In this discrete formulation of a gauge theory, belonging to a compact 
group G, the gauge potential is replaced by a map of bonds (ordered pairs 
b of nearest neighbor lattice points) into G: b ^ gi, & G, with the property 
that Qx^y = Qyl for any bond b = {x, y) . In order to formulate the analog 
of the Yang-Mills action we need the notion of a plaquette. This is a closed 
curve of four bonds P = {x, y), {y, z), {z, u){u, x). To each of these we assign 
the group element gap — gx,y9y,zgz,ugu,x, a-nd for a given unitary character x 
of the group G the real number 



l{x{gdP)+X%gap))- (39) 



Note that Sp is independent of the orientation of P. More generally, we 
can assign to each closed path C on the lattice the product gc of its g^s, 
along one of the orientations, and then take its character x(C). This number 
is also independent of which site we begin the loop with. The action of a 
configuration {g} on a finite region A of the lattice is taken to be 

SA{g) = -T.Sp. (40) 

PCA 

The expectation value of an 'observable' A{{g}) in A is defined by 

(^)a = ^J exp(-5A)A diiA, Za = J e-^^ diiA, (41) 

where dii\ is the product of the normalized Haar measures for all bonds in 
A. 

A local gauge transformation is a map x '-/^ & G. Under such a 
transformation a gauge field configuration transforms as 

gx,y ^ lxgx,yly^ {b = {x, y)). (42) 

Clearly, Sp, the action and the measures d^\ are invariant under local gauge 
tr ansf ormat ions . 

Using the Baker-Hausdorff formula, it is easy to see that the naive formal 
continuum hmit leads to the Yang-Mills theory. There are good reasons to 
believe that the continuum limit of a pure lattice gauge theory provides a 
consistent interacting Euchdean field theory. 

Next, we add in a gauge invariant manner Higgs fields to the lattice model. 
To do this for such a field 4>{x), we have first to introduce a discrete form of 
the covariant derivative. If this field transforms with respect to G according 
to the unitary representation U (g), then the following definition has the right 
formal continuum limit (we use the notation introduced in Appendix B) 

Dn(p{x) = -\u {g:c,x+aei,)(t>{x + oc^) - 0(a;)] . (43) 
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The adjoint operator is 

Dl4>{x) = -[u {gx,x-ae^)4>{x - flC^) - . (44) 

The corresponding gauge invariant Laplace operator (generahzing (57) in 
Appendix B) is 

= -DlD,. (45) 

The discrete Higgs part of the action is then in D dimensions (dropping the 
A-dependence) 

S^[{g}, 0] = Y: I Ux,y<Piy) - m ? ^ \/(0(x)), (46) 

{x,y) X 

where U^^y := U{g^^y). 

When fermionic fields are discretized in the most obvious way, the prob- 
lem of 'fermion doubling' shows up. Early remedies to cure this involved 
an explicit breaking of chiral symmetry of QCD (see, e.g.. Chap. 4 of jH]). 
Fortunately, there has been recently important progress in solving this long- 
standing problem. Lattice formulations that preserve chiral symmetry have 
been found, but in practice numerical simulations become much harder. For 
a very readable recent review on these aspects, see jH]. 



Continuum limit 

What we are really interested in is what the theory describes in the continuum 
limit. The lattice is only used as an ultraviolet gauge invariant regularization, 
valid outside the domain of perturbation theory. 

Consider, for concreteness, a pure Yang-Mills lattice theory with a single 
dimensionless (bare) coupling constant g. The only dimensionful parameter 
in the theory is the lattice spacing a. Therefore, a physical quantity with the 
dimension of a mass, say, such as a particle mass rriphys, must be of the form 

rriphys = -f{g)- 
a 

The theory allows us to compute /((?), and thus ratios of masses. Clearly, in 
taking the continuum limit a — *• we have to tune g a,s a, function of a such 
that aruphys 0. Since the inverse, ^ = {amphys)~^ , is a correlation length, 
we have to approach - in the language of statistical mechanics - a critical 
point of the theory. The tuning of g{a) is controlled by a renormalization 
group equation of the form 

= -PhAT^g)- 
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In lattice perturbation theory (3lat{9) can be computed as a power series in 
g for small g: 

m = -f^og' - f3ig' + .... 

For SU{N) one finds 

N 11 
Po = TT^^' etc. 

From this we find by integration 

2_ 1 

^ /5oln(a2Ai^^)^-' 
where A^at is an integration constant. Inverting this one sees that 

/ 1 

aruphys = const exp(^-— — 



showing the non-perturbative character of masses in the theory. Note that 
the bare coupling constant g{a) vanishes in the continuum limit. 

The dimensionful lattice A-parameter, A^at, provides a scale which sur- 
vives the continuum limit. This is remarkable, because classically the theory 
has no scale. Through the process of renormalization we have introduced a 
mass scale into the quantized theory, a mechanism that is often called di- 
mensional transmutation. This renormalization procedure is in the lattice 
formulation no more mysterious than adjusting the temperature toward the 
Curie point of a ferromagnet. At the same time these considerations shows 
that the continuum limit has to be taken with care. They also show an in- 
teresting point of contact with critical behavior, discussed in other lectures 
at this school. 



Elitzur's theorem for gauge theories 

For systems with global symmetries the phenomenon of spontaneous symme- 
try breaking (SSB) is well-known. This is accompanied by a non-vanishing 
'magnetization'. At first sight one expects something similar for gauge the- 
ories. However, Elitzur |B] has shown that local quantities, like the bond 
variables gt or a Higgs field, which are not gauge invariant, have always 
vanishing mean values. This is quite easy to prove, and is also physically 
understandable. It is instructive to show this first for a lattice gauge theory 
with gauge group Z2, because this shows the contrast to the Ising model. In 
Appendix C we give a general proof of Elitzur's theorem. 

The bond variables will be denoted by ab{= ±1. The action for a finite 
region A of the lattice is taken to be 

SAi{a}) = -J2^0P-hT.^b, (47) 

PcA 6cA 
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where h is an external 'field'. We could also include a Higgs field trans- 
forming according to i—> £j(/)j, £j = ±1. 

Before continuing with this model, let us recall the SSB for the Ising 
model with its global symmetry Z2, consisting of the identity and the reflec- 
tion (Tj I— s> — (Tj for the Ising spins cXi of all lattice sites i. Above a critical 
temperature there is only one infinite-volume equilibrium (Gibbs) state. 
However, for T < Tc each translation invariant equilibrium state (proba- 
bility measure) fj,^ is a convex linear combination of two different extremal 
states This means that fi^ describes a mixture of two pure phases. The 
probability measures are weak limits of Gibbs states on finite regions 
A C with ± boundary conditions outside A. Since they are different, 
they are not invariant under the symmetry group Z2 of the interaction; the 
symmetry is spontaneously broken for these pure phases. Correspondingly, 
the spontaneous magnetizations 

= (a,)^, (48) 

do not vanish for j3 > j3c. 

In sharp contrast to this situation, the mean value of (ub) does not signal 
a symmetry breaking for the Z2 lattice gauge model: 

Theorem (Elitzur). For the expectation value {ab)\{h) we have 

lim{ab) — uniformly in A and (5. (49) 

In particular, the thermodynamic limit of {o'b)\{h) vanishes for h I (no 
spontaneous 'magnetization'). 



Before giving the simple proof, we remark that when a Higgs field (f)x is 
added, a similar proof shows that 

= (50) 

(Exercise) . This does, however, not exclude a Higgs phase with an exponen- 
tial fall off of the truncated correlation function {(TQp^aQp^)c for the gauge 
fields (mass generation), but (49) and (50) show that this is not signaled by 
local gauge variant observables. This fact is well-known to people working in 
lattice gauge theory, h\it is largely ignored outside this community. As order 
parameters one may try to choose local gauge invariant quantities, such as 
the norm of a Higgs field. 
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Proof of Elitzur's theorem. We can choose in (46) the bond variable cxoi 
for the bond b — (0, 1), and estimate in 

, > Y.{a}0-oiexp{/3Y.pO-dP + /3hY.b(^b) 
exp(/3 Ep (^dP + Ph J2b (^b) 

the numerator and the denominator D separately. 

Consider a gauge transformation (T.-^- = eiaijSj, with £j = 1 for i 7^ 
and replace aij in N and D by Eia'^jCj, dropping afterwards the prime of 
(7y {aij I— > SiUijej). This can be done for = il- D is equal to half the 
sum: 

2D 

D ^ - ^eo^e-x.^{j5^aQp + (5hY^ Vb) exp {eo(5h ^ a^j) , 

{cr} P b j=l 

where the prime of the sum means, that the bonds (0, j) must be excluded. 
Clearly, for /i > 

L> > ^ exp(/3 ^ agp + ^ V,) le''^'^\ 

{a} P b ^ 

Similarly, we have for the numerator 

^ 2D 

MP 6 ^ j=l 

and thus 

|iV| < 5]cxp(/3^a9P + /3/i^V(,)sinh(2D/3/i). 

{a} P 6 

This gives the estimate 

\{(Toi)A{h)\ < 2e^^'^^smh{2D(5h) ^0 for /i ^ 0, uniformly in A and /3 . 

(52) 

What is the physical reason for this different behavior of models with local 
and global symmetries? Consider once more the Ising model in the absence 
of an external field. At low temperatures, the two regions in configuration 
space with opposite magnetizations and — cr^ can only be connected by 
a path involving the creation of an infinite interface which costs an infinite 
amount of energy. Therefore, the process cannot occur spontaneously. Al- 
ternatively, one can make use of a small external field which is switched off 
after the thermodynamic limit is taken. In this procedure, one of the two 
fundamental configurations {a.^ = +1} and {a^ = —1} is energetically elim- 
inated and a non-vanishing magnetization is left. On the other hand, in a 
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local gauge theory one can perform symmetry transformations which act only 
non-trivially on a finite set of basic variables on which a local 'observable' 
depends. The system behaves therefore similarly as quantum mechanical 
systems with finite numbers of degrees of freedom. This aspect is manifest 
in the proof given above, and becomes perhaps even more transparent in the 
general proof of Elitzur's theorem given in Appendix C. Sometimes it is said 
(see, e.g., p. 810 of Ref. jlj) that the global symmetry associated with the 
gauge group is spontaneously broken. The precise meaning of this statement 
is, however, unclear to me. 

Gauge invariant effective potential 

From lattice simulations it is known that the expectation value of |0p is 
suited to characterize the Higgs phase. It is, therefore, natural to couple 
the source in the partition sum to this quantity, and consider the effective 
potential belonging to 



where djj^A is now the product of the normalized Haar measures and the 
Lebesque measure [d(f)]\. We emphasize that this effective potential should be 
regarded as a function of since this is the conjugate variable, and only as 
a function of this variable it has to be convex. (This remark is, I believe, also 
relevant for the Ginzburg-Landau free energy density of superconductivity, 
which is only convex as a function of the square of the absolute magnitude 
of the order parameter. Drawing this free energy as a 'Mexican hat' of the 
complex order parameter is somewhat misleading.) 

For gauge models with the sum of the actions (40) and (46) the effective 
potential can be computed numerically or in perturbation theory. The latter 
has been developed on the lattice in close analogy to the continuum gauge 
theory [?]. We shall elaborate on this in Sect. 3.1. 

Appendices to Chapter 2 

Appendix A. The Legendre-Fenchel transform 

In thermodynamics and statistical mechanics it is important to generalize the 
usual Legendre transformation to functions which are not everywhere differ- 
entiable. This generalization in a finite number of dimensions was introduced 
and developed by Fenchel in 1948. 




(53) 
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The Legendre-Fenchel transform f* of a real valued-function / on i?" is 
defined as 

f*{x*) = snp{{x,x*) - fix)}, x*eK'. (54) 

To be precise we have to assume that / majorizes at least one convex function. 
(This definition has an immediate generalization to dual pairs of infinite 
dimensional spaces.) The convex hull of / is then defined as the largest 
convex function majorized by / and is denoted by cl{conv /). In this situation 
/* is convex and the following statements hold: 

/* = {diconv /))*, /** = diconv /). (55) 

For a proof and further information we refer to p2]. For a useful treat- 
ment, see also [13] or [Hj. 

If / is convex and differentiable, this generalized Legendre transformation 
reduces to the standard one (show this). In one dimension there is a simple 
geometrical construction of /* (exercise). 

Appendix B. Computation of a functional 
determinant 

In this Appendix we compute the functional determinant of the operator 

A = -A + m^ + P"{^) (56) 

that appeared in (23). It may be instructive to start from the lattice regular- 
ization and then perform a (formal) continuum limit. The discrete Laplacian 
is naturally defined by 

-A = J2did^, (57) 

a 

where 

idj)ix) = ^[fix + ae„) - fix)], idif)ix) = ^[/(x - ae^) - fix)]. (58) 

Here are the unit vectors in the lattice directions, and a is the lattice 
constant. The Laplacian becomes diagonal in the Fourier representation. 
On a finite sublattice with Na lattice points in the direction Cq,, the Fourier 
transform of a function fix) is normalized as 

fik) = Y: fix) e-^'-\ (59) 
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The inverse transformation is given by the following sum over the discrete 
first Brillouin zone A ^ {k : ka = (7r/a)(2na/iV„), - Nc,/2 < n„ < Ar„/2}: 



H^l fceA 

One easily finds for the Laplacian in Fourier space: 
-A/^ ^2(l-cosM;(,). 



(60) 



(61) 



The Fourier representation of the operator A is therefore in D dimensions 



D 



a=l 



Hence, 



2 

In det A = J2 1^ [-^^ + — ~ C'^s ^«^) 



(62) 



(63) 



where = + P"($). In the thermodynamic limit Ni = N2 = —: 
A?" — > 00, we obtain, using the rule 



fcGA " -"^ fe4 

with BZ = [-^, 

IndetA /■ c?^^ v^, 2 

= / 7 — > In + — > (1 - cosA;aC 
(A^a 3 Jbz (27r 3 4^ L 



(64) 



fc4 



Here, A;4 = (27r//5)n, -/5/2a < n < /?/2a. 

This quantity diverges, of course, in the continuum limit a — > 0. We now 
show that the temperature dependent part of the regularized free energy 
density 

11 

converges to (24). Note that this formula reduces for P = to the correct 
free energy density of a free scalar field with mass m. 

For a — > we can replace the argument of the logarithm in (64) by 

n^{k) + kl Thus 



freg 



2/3 Jbz (2%) 



n=-l3/2a 
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(65) 



From this we subtract its value for /3 — > oo. The sum appearing in (35) can 
be worked out with the ^-function technique. We use 



where 



/3/2a 



2^2. 



n=-l3/2a 



ds 



\s=0 ' 



1 



l3/2a 





For a — > 0, (^^^ converges to 



(66) 



(67) 



+ 



/32 



(68) 



This is the ^-function of the 1-dimensional differential operator — ((i/rfr)^+n^ 
on the interval [0, /3] with periodic boundary conditions. Subtracting from 

(^'"^■''(.s) its vahic for /3 — >■ oo amounts for its derivative at ,s = to the 
subtraction of a term linear in and thus of a /9-independent term for /'"'^s. 
Hence we need 



13 a^O 



reg 



ds 



d 

ds 



for s = 0. 

In the footnote below we show that 



C'(0) = -/3Q-21n(l-e-^^). 



(69) 



This proves our claim. Note that the first term in the last equation gives the 
zero-point energy. 



Footnote: Calculation of the C-function 

We recall that the ^-function of a self-adjoint positive operator Q with a pure 
point spectrum and the spectral representation 

n 

where P„ are the projections on the eigenspaces, is defined as 

Cq{s) = E fl' dim P„. (70) 

s 
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From det Q = exp(Tr InQ), we get 



In det Q = ^ In A„ = - 



dCq 



ds 



(71) 



s=0 



Note also 



1 roo 



1 f°° 



r{s) 



roo , . 

dtf-^Tr{e-*^). (72) 



For the operator Q — —{d/drY + o;^ on the interval [0, /3] with periodic 
boundary conditions, we obtain 



n=-oo [/32a;2 + 4^2^: 



r(.) 



/ daa^-^ V exp -a(/3V + 47rV 



n=— oo 
\2Z 



= ^-^^ + ^i:^(-l)'E r da a^^^-'e-'-'-'^ 



n=l 







where Cr{s) is the Riemann ^-function 

oo -1 



n=l 



Since Cii(O) = -|, C^i(O) = -|ln(27r) we have 

C^(0) = -21n(/?^) + 2f Ci^(20- 
The last factor can be expressed in terms of the Bernoulli numbers B21: 



CK(20^(-l)«|g 



B. 



21- 



(73) 



Therefore, 



C^(0) = -21n(/3u;)-E^Ti^2^ 



Here we use the well-known relation 
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Integrating leads to 

Using this we obtain indeed (39): 
(^(0) = -2 \vl{(5uj) + 2 In ^ - 2 In sinh ^ = -(5uj - 2 ln(l - e'^'^) . (74) 

Appendix C. General proof of Elitzur's 
theorem 

It is quite easy to show generally that in a lattice gauge theory the expectation 
value of any function, which depends only on the basic variables in a finite 
region (local observable) and is non invariant under the gauge group, has 
a vanishing expectation value, when a source term is set to zero after the 
thermodynamic limit is taken. 

In what follows, the basic dynamical fields (gauge fields gx^y, Higgs fields 
(t)x, etc.) will collectively be denoted by (/?. As in (44), we subtract from the 
gauge invariant action S[(p\ a coupling J ■ to a. (small) external source field 
J. The expectation value of an 'observable' A[(p] in a region A is, as always, 

{A)k,j = ^ / exp(-5A + J-^)A dfXA, Za,j = J e-^^+'^-^ d/XA, (75) 

where d/iA is the standard gauge invariant measure (involving the Haar mea- 
sure for gx,y, etc.). We denote the action of a gauge transformation g of the 
basic fields by (/? i-^ T{g)(f. Non-invariancc of a local A[lp], depending only 
on field configurations belonging to a finite subset of the lattice, means that 

J A[T{g)ip] diiA = 0, (76) 

because the integral on the left gives the invariant component of A[ip\. 

We split the field configurations into two complementary sets {(/?'} and 
where the first set denotes those on which A depends. Associated 
to this we consider the subgroup of gauge transformations which keeps the 
elements of the second set fixed. If {dg)x denotes the product of the Haar 
measures corresponding to this subgroup, we have 

{A)a,j = ^ UM-SM + J' ■ T{g)ip' + J" ■ ip")A[T{g)ip'] di,A{ip){dg)A. 
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Because of (72) we can subtract 1 from the factor exp[J' ■ T{g)ip'] in the 
integrand, hence the integral vanishes for J — > 0, for 'reasonable' observables 
A, uniformly in A and (3. Since the partition sum remains uniformly positive 
in this limit, we arrive at 

(A) = lim^lim^(A)A,, = 0. (77) 

Intuitively one can perhaps say, that a gauge theory cannot be sponta- 
neously broken, because local symmetry transformations involve only "local 
changes", and hence the transformed states are not disjoint, similarly as 
for systems with a finite number of degrees of freedom. Unfortunately, the 
commonly used terminology in most texts does not describe the situation 
adequately. 

3 Applications to cosmological phase 
transitions 

We now turn to relevant applications of the general theory to cosmological 
phase transitions. 1 shall begin with a relatively detailed discussion of the 
electroweak phase transition, on which much effort gas been invested in recent 
years. 

3.1 Thermodynamics of the electroweak phase transi- 
tion 

One of the main motivations for the large amount of work was the suggestion 
that the observed baryon asymmetry in the Universe was eventually deter- 
mined at this transition. Indeed, all three Sakharov conditions for baryo- 
genesis are, in principle, fulfilled. One of them, namely the departure from 
equilibrium, depends on the nature of the phase transition. 

Out of equilibrium processes require a first order phase transition. We 
shall see that this would be reahzed for Higgs masses much below the current 
experimental lower limit. With increasing Higgs mass the phase transition 
weakens rapidly and presumably turns into a smooth crossover. For realis- 
tic Higgs masses it becomes too weak for baryogenesis. In supersymmetric 
extensions of the Standard Model there remains a small parameter range for 
which the transition is sufficiently strong. 
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3.1.1 The 5f/(2)-Higgs model 



Reliable results for the electroweak phase transition can only come from non- 
perturbative methods, that is from lattice simulations, because the Higgs 
field is nearly massless at the transition which leads to serious infrared di- 
vergencies in the perturbation expansion. A complete lattice simulation of 
the Standard Model is not feasible due to the presence of chiral fermions. 
The main (infrared) problems are, however, only connected to the bosonic 
sector. For this reason people usually studied simulations of the S'f/(2)-Higgs 
model and used perturbation calculations to include the U{1) gauge group 
and fermions. 

In what follows, I shall describe some results of such studies. For practical 
reasons (avoidance of huge lattices), people used in recent years asymmetric 
lattices (different spacings in temporal and spatial directions; see e.g. ^IS^). 
For the sake of presentation, I shall ignore the complications connected with 
this. 



The lattice action 

In writing down the action we adapt the previously used notation to what 
has become standard.^ The bond variables of the gauge group SU{2) are now 
denoted by Uxfj, for the bond {x,x + fi), where fi stands for unit vector e^. 
The Higgs doublet can be described by a 2 x 2 complex matrix ipx, satisfying 
the pseudo-hermiticity condition 



t T -1 



1 
-1 



Under a gauge transformation x ^ Vx the Higgs field transforms as ip^ ^ 
Vxfx- Equivalently, the Higgs doublet ca be described by four real scalar 
fields 05, related to ip by 

f = 4>oh + i4>kn, 

where the are the Pauli matrices. 

In the formulae below we set the lattice spacing a equal to unity. The 
action is parametrized as follows: 

S[U,^]=(3Y.i^-l-TrUp) 
p ^ 



2 



(1 rl ]2 ^ 



^Slight remaining differences to some of the quoted papers amount to the replacement 

u —f r/t — r7~i 
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Here, the parameter (3 (not to be confused with l/kT) is related to the 
gauge couphng constant g by /3 = 4:/g'^- The last term is the hopping part 
of the gauge invariant kinetic contribution of the scalar field. The hopping 
parameter k is related to the bare scalar mass mo by itLq = (1 — 8k — 2A)/k 
(show this). As mentioned earlier, the gauge invariant length of the Higgs 
field, := ^Tr{(pl(px), plays an important role. It turns out that this 
quantity has large autocorrelations and can be used to characterize the Higgs 
phase (as in the standard semiclassical treatment of the Higgs mechanism). 

Gauge inveiriant effective potential 

An important and useful object for investigating phase transitions is the 
gauge invariant effective potential (53). For the present model we have to 
start from the partition sum 

= y e-^^+"^S-eAP' dp A, (79) 
with the action (78) which we rewrite, using the polar decomposition 

<Px = PxOix: ax e SU{2), (80) 

as 

S[U,^]=Pj2{l-lTrUp) 
p ^ 

4 

+ E{(P-)' + KpI - 1)' - «EP-Pa:+ATr(a,)^t/,^a,+^)}. (81) 

X fl=l 

The expression in the last trace is a gauge invariant link variable, which we 
call 14^. Because of the gauge invariance of the gauge action Sg[U], we find 
for the total action 

S[V, p] = S,[V] + J2{pI + KpI -^f-'^T. PxPx+f.Tr{V,,)}. (82) 

X n=l 

Next, we have to discuss the integration measure. Prom the link variables 
we get the product of the Haar measures dV^^. For the Higgs field we 
have from every lattice point the product of the Lebesque measures for the 
real fields, 

d^x = ]Jd(ps{x) = pldpxda^^, 
s 

where da^ is the natural measure of the three sphere in 0-space. Using (80) 
this can be written as 

dipx = pldpx^Ti'^dax, 



28 



where da^ is again the Haar measure. Since the action does not depend on 
CKj; we can trivially integrate over the da^ and obtain for d/i in (79), up to 
an irrelevant normalization, 

dii = l[dV,^l[pldp,. (83) 

X,IJ, X 

Eqs. (79), (82) and (83) determine the thermodynamics of the SU{2) Higgs 
model. Note that - without any gauge fixing - three scalars have disappeared. 
(Repeat this in the continuum theory.) The action (82) and the measure (83) 
have a global 5" C/(2) invariance: px ^ Px, V^fj, i— > UVx^U'^, U e SU{2) {weak 
isospin). This symmetry is reflected in the mass spectrum. 

The phase structure of the model will be discussed shortly. Let us consider 
a first order phase transition between the 'symmetric' (s) and the 'broken' 
(b) phases. Along the coexistence curve J{T) the free energy densities are 
the same: 

Ws{T, J{T))^w,{T, J{T)). 

Differentiating with respect to T gives (note that w{T^ J) is the negative of 
the Hclmholtz free energy density) the following Clausius- Clapeyron equation 
between the latent heat AQ and the jump Ap^ in the order parameter 

AQ = Ap^T^. (84) 



Phase structure 

Numerical studies established for a fixed A = 0(1) the following picture (see 
Fig. 4) of the phase diagram in the 3-dimensional parameter space (3, k, T, 
where T is the temperature in lattice units (inverse of the time extent of the 
lattice): 

For K, = the model reduces to a pure SU (2) gauge theory. As indicated 
in Fig. 4, there is for this limiting model for any T > a deconfinement phase 
transition at a critical value /3c which increases with decreasing temperature. 
This phase transition also exists when k is switched on, but for large values of 
the hopping parameter the transition may change to an analytical crossover. 
In the limit (3 = oo the gauge fields become trivial (pure gauge). Setting them 
to unity, we get a (p^ theory for an SU (2) Higgs field, having a global SU (2)^ x 
SU{2)ji symmetry. For a fixed positive temperature this model exhibits a 
phase transition for a critical value Kc, above which the symmetry is broken 
down to 0(3). Simulations showed always a second order phase transition. 
With varying temperature a transition line is found. This develops into a 
surface for finite values of /3 where the electroweak phase transition takes 
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Figure 4: Phase diagram for the SU{2) Higgs model. For a fixed quartic 
couphng A = 0(1) there are three regions in the /?, k, T parameter space: 
confinement region for f3 < f3c, <^ 1; Higgs region for ^ 1, k, > Kc] 
deconfinement region for (3 > jSc, k, < Re- 



place. When /3 becomes small, one expects a region in parameter space 
where no phase transition occurs. So far A was kept fixed. When both A and 
the temperature become large enough, the electroweak transition may turn 
into a crossover phenomenon. 

For T = the phase structure is sketched in Fig. 5 in the parameter 
space [3, k, A. For small values of k the system is in the confinement phase 
(at K = the model reduces to a pure SU{2) gauge theory). Along the 
surface S there is a first order phase transition to the Higgs phase (with one 
spin-0 Higgs boson and three massive gauge bosons), except at the boundary 
(5 = oo {g = Q). The boundary dT^i of is a second order transition line 
as for T > 0. For small {3 (large coupling) and relatively large values of the 
hopping parameter k there is an analytic connection between the two regions 
beyond the boundary piece dTi2- (This has been established rigorously in 
|16j.) The distinction between confinement and Higgs phases has then no 
meaning anymore. 

If one uses the actual value of the weak coupling constant, it turns out 
that the finite temperature phase transition always takes place between the 
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X= CO 



p P=oo 

Figure 5: Phase structure for zero temperature (adapted from [7], Sect. 6.1). 

Higgs and the deconfinement regions. 



One-loop lattice calculation of the effective potential 

The effective gauge invariant potential of the SU{2) Higgs model has been 
determined numerically, along with other quantities jl2j. Here, we want 
to derive the one-loop approximation of this important object. This is an 
instructive example of lattice perturbation theory. 

Let us first consider the 'broken' phase. The exponent in (79) has for 
= 1 a stationary point for the value p given by 



, 8k + J - 1 



i5) 



provided J > 1 — 8/t — 2A. In the tree approximation (Laplace approximation 
of the infinite volume integral) the effective potential is 



Vtreei^^) = (1 " S/t)^^ + A(<1>2 - 1)^. 



(86) 



(This is convex as a function of $^ !) For the one-loop correction we set 
Px = P + o'x, = exp(A^(x)) and consider A^{x) as small. The quadratic 
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piece in of the exponent in (79) is just the action of a free scalar lattice 
field a with effective mass 



vr? = -\p^ = mH,Qp^/v'^, (87) 

where v is the value of |$| at the minimum of Vtree(*^*^), and mH,o is the 
zeroth order Higgs mass. According to Sect. 2.2 this gives the following 
contribution to the one-loop correction 

where k"^ := X]^2(l — cosA;^). For T > the integral over has to be 
replaced by T times the sum over k4, with k^ = (27rT)n, n G Z. 

For the contribution of the gauge field we set V^^ = exp(y4^(a;)) and 
consider A^(x) to be small. Then the gauge part of the action reduces in 
quadratic approximation to the Maxwell action of the three scalar fields A"^, 
defined by = These receive, as in the continuum theory, a mass 

term from the coupling to the scalar field (last term in (82)) with effective 
mass 

1 

2 

where mwfi is the zeroth order W-mass. 

Next the measure (83) has to be rewritten in the appropriate approxi- 
mation in terms of Afj_{x) and Ux- It is a nice exercise in group theory to 
express the normalized Haar measure dU of SU (2) in terms of the canonical 
coordinates (of the first kind) Xa, defined hy U — exp(iXaTa/2). The result 
is 

Hence, 

1 sin^(,|A,|/2) 

with \A^\^ = We write the first factor in (83) as 

J{dV,^ = e-'-J{dAl, 



2 -^2-2 2-2/2 foc\\ 

^g=n^^9P = m^fiP h > (89) 



with 
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(An additive numerical constant was dropped in the last expression.) 

In sufficient approximation the measure (83) becomes, up to a normal- 
ization, 

d^^ = exp (- Y: ^AI.A;] n dAl{x)l[da, . (92) 

The exponential in this expression is a mass term with mass g^/Q. This acts 
like a mass-counter term and diverges in the limit a ^ 0. We encounter here 
an example of a general phenomenon: The lattice regularization provides its 
own counter terms which ensure renormalizability of the theory. 

The contribution of the free massive vector fields An{x) with the effective 
mass rn^g + gives the following contribution to the one- loop potential 
(exercise) 

= 9 L ^ Hk' + m'w^Vv' + //6). (93) 

The reason for the factor 9 is obvious: three massive vector fields with three 
degrees of freedom. The additive term g'^/ 6 in the argument of the logarithm 
would not appear in a continuum calculation. 

We consider now the 'symmetric' phase. There is a stationary point at 
p = 0, for which the tree level potential vanishes. The exponent in (79) 
is now (apart from a sign) the action of a free scalar lattice field px with 
effective J-dependent mass — ml — J/k, plus that of three massless 
vector fields A^{x) which get an effective mass g'^/Q from the link integration 
measure in (83). We have to keep in mind that is a radial variable and 
that the measure (83) contains the weight factors p^. In the steepest decent 
approximation we obtain for the one-loop free energy density, apart from 
J-independent contributions due to the gauge fields, 

+ (94) 

The one- loop effective potential is the Legendre transform of w{J) ($^ is the 
conjugate variable). This has to be done numerically. 

At the critical temperature the minima of the two effective potentials are 
degenerate. The critical temperature in units of the Higgs mass as a function 
of the Higgs mass is shown in the upper part of Fig. 6, both for numerical 
simulations and ])erturl)ation calculations.^ 

"^We do not discuss here the so-called effective field theory method, which makes use of 
the compactness of the time interval for positive temperatures. This suggests, in complete 
analogy to the dimensional reduction in Kaluza-Klein theories, the construction of a 3Z) 



33 



1.5 



J I I I I I I L. 



J I L 



20 40 60 



80 



T — r 



T — r 



T — r 




20 40 60 80 
M^LGeV] 

Figure 6: {top panel): Critical temperature as a function of the Higgs mass 
tuh- The full squares and triangles are from 4-dimensional simulations with 
different lattice extensions in the imaginary time direction. The full circles 
are results from lattice perturbation theory, extrapolated to the continuum. 
The solid lines are from continuum 2-loop calculations fH]- (From Ref. [TH].) 
[bottom panel): Jump of the order parameter at the critical temperature. 
The triangles and squares are from the same simulations as in the top panel, 
and the solid line is again the result of a continuum 2-loop calculation. The 
starred symbols are from simulations of the so-called 3-dimensional 'reduced 
model' (see footnote 4). 
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Figure 7: Latent heat as a function of the Higgs mass. The symbols have 
the same meaning as in the previous figures. 



The jump A in the gauge invariant quantity pi, vt '■= (A(p^))^/^, is a 
measure of the strength of the first order phase transition. Its dependence 
on the Higgs mass is shown in the lower part of Fig. 6. A similar plot for 
the latent heat is given in Fig. 7. More recently it has been shown ^S] 
that the first order transition ends at mn = 66.5 ± 1.4 GeV. Above this 
endpoint only a rapid cross-over can be seen. This number increases for the 
Standard Model to 72.1 ± 1.4 GeV. In the framework of the effective field 
theory method (footnote 4) this result has been obtained in [21], |22| and 
122 (for a short review, see ^^)- 

3.2 Baryogenesis and electroweak phase transition 

The previous results, in particular Fig. 6, are crucial for answering the ques- 
tion whether the cosmological baryon asymmetry might have been generated 
in the electroweak phase transition. We shall see that this was not the 
case. Nevertheless, fermion-number violating electroweak interactions are of 

effective action containing the zero modes of bosonic fields. For an introduction to this 
method ( which can be regarded as a generalization of the high temperature expansion), 
as well as for extensive references, see PT)] . 
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interest, since they changed the baryon and lepton asymmetries that were 
generated during earher epochs, expected in the context of unified extensions 
of electroweak and strong interactions (GUTs). For reviews, see j^Hj and [26] . 



Anomalous fermion number violations 

Let me first recall why baryon number (B) and lepton number (L) are not 
conserved in the Standard Model, although they are obviously conserved in 
perturbation theory. On the non-perturbative level this due to the chiral 
nature of the weak interactions. Indeed, this implies that the divergences 
of the baryon and lepton currents are both anomalous (quantum mechani- 
cally non-conserved). It turns out that this violation is such that B — L is 
conserved. 

Consider, for illustration, in the SU{2) Higgs model a massless left- 
handed fermion doublet ipL and the global U{1) fermion current 

Classically, this is conserved, and expresses the conservation of the number 
of fermions corresponding to the field ipi- In quantum theory, the famous 
triangle anomaly of Adler, Bell and Jackiw leads to the modification^ 



^A'^'' = -J^HF'.uFn, (95) 



where is the dual field strength of the SU{2) gauge field F^^, belonging 
to the gauge potential Af^{x) = (see, e.g.. Chap. 19 of [2H1 or 

)• The right-hand side of this equation, the famous characteristic second 
Chern density, has an important topological meaning. Its integral (Chern 
number) is equal to a certain winding number, as shown below. Locally, 
it is a divergence of a current. It is more elegant to write this in terms of 
differential forms. The Chern 4-form belonging to the right-hand side of (95) 

C2{A) = ^tr{F A F) {F = ^F^.dxf' A dx") (96) 
is (modulo global questions) an exact differential 



1 



87r2' 



2 



C2 = dQ{A), Q{A) = —Tr AA{dA + -AAA) . (97) 



3 



^Thanks to the work of M. Liischer (30], this can now also be derived on the lattice 
in complete analogy to the Fujikawa method for the continuum theory. For review and 
references, see [O]- 
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Q{A) is the Chern-Simons form {A = A^dx^). This 3-form is not gauge 
invariant. If we integrate (95), i.e., 

d^J = C2iA) (98) 

over a region V of space-time between two time shces. Stokes' theorem gives 

[ di^J = I i^J - I = ANcs, 

JV Jt=t2 Jt=tl 

where Ncs is the difference of the Chern-Simons numbers 

ANcslA] = [ Q{A) - [ Q{A). (99) 

Jt=t2 Jt=tl 

This difference is gauge invariant since it is equal to the integral of C2{A) 
over T>. We see that the change of the fermion number ANp associated to 
is equal to ANcs- 

ANf = ANcs[A]. (100) 

Before further elaborating on this, we consider the left-handed doublets of 
the three families of leptons and quarks of the Standard Model. The previous 
equation implies the same fermion-number change for each species. For the 
leptons the changes of lepton number are 

ANe = AN^ = AN, = ANcs[A], (101) 

and the changes for the quarks imply the following change of the baryon 
number 

= ^ • 3 • 3 ANcs[A] (102) 

(| from the baryon number of quarks, 3 • 3 from the three colors and number 
of generations). This also shows that B — L is conserved. 

The non-perturbativc effects we arc discussing involve strong gauge fields, 
and it is therefore reasonable to treat them classically. So let us now consider 
classical field configurations and the changes of their Chern-Simons numbers 
in time. In particular, we are interested in field configurations which interpo- 
late between different classical vacua with pure gauges = C/~^5^C/, U{x) 
an S'f/(2)-valued function. The Chern-Simons number of such a vacuum 
gauge field is in general different from zero. A short calculation gives 

Q(A) = --L^Tr(u-^dU A U-^dU A U-'^dU). (103) 

We assume that for fixed time U{x) converges at infinity to I2 in such a 
manner that we can regard it as a function on the compactified space S^. 
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Since SU{2) is topologically also a three-sphere, we then have a map from 
5''^ to S^. We now show that the winding number (degree) of this map is just 
the Chern-Simons number of the pure gauge. The easiest way to see this is 
to note first that U~^dU is the pull-back by U of the Maurer-Cartan form 6 
on SU{2y. 

U-^dU = U*{9). 
Therefore, Eq.(103) can be written as 

Q(U-^dU) = --^U*\Tr{e AO A9)]. (104) 
247r^ L J 

A general result on winding numbers tells us that the one for U : 

for any volume form u on 5*^ (see, e.g.. Sect. 7.5 of [3l|), in particular for 
cu = Tr{9A9A9). Hence, 

Ncs[U-'dU]=n{U){--^) f uj. 

247r2 Jsui2) 

But (— l/247r^)c<j is the three-form corresponding to the normalized Haar 
measure, whence 

NcsP'^dU] = n{U). (105) 

This integer characterizes topologically different vacua. 

As a result of these considerations we obtain for a field configuration 
which interpolates between vacua at t = — oo and t = +oo 

ANcs[U-'dU] = niU)\,^^ - n(f/)|,^_^ . (106) 

This difference is also equal to the second Chern number 

C2 = I C2(A), 

which is, of course, gauge invariant. 

The winding number n{U) also shows up in the vacuum configuration of 
the Higgs field. It is straightforward to verify the following identity (we omit 
the wedge symbols): 

Tr[{(f)U(f)f]-3Tr[A{dA+'^A'^)] = Tr[{cf>^D(t>f]-3Tr{^]FD(t>)+3d{Ad(f)(f)^). 

(107) 
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This implies that the left-hand side is gauge invariant, up to an exact differ- 
ential. Apart from this exact differential the right-hand side vanishes for a 
classical vacuum, where Dcf) = 0. Integrating the last equation we see that 
the winding number of the Higgs vacuum field 

= ^ / Tr[(0trf0)3] (108) 

is equal to the Chern-Simons number, and for a vacuum-vacuum transition 
we obtain the gauge invariant relation 

ANcs = ANh. (109) 

For the smallest non- vanishing change (102) gives 

AB = -AL = 3, (110) 

the quantum numbers of 12 left-handed fermions. States such as YluLdhdhi^L-, 
where 11 denotes the product over all generations, can be created out of the 
vacuum. (Altogether there are 12 such combinations.) 

Vacua with different topological numbers n are separated by a potential 
barrier, because a transition has to involve non-vacuum fields. At zero tem- 
perature and 'low' energies such a transition is a tunnelling process that can 
be estimated by using (constraint) instantons. For a pure Yang-Mills the- 
ory an instanton minimizes the Euclidean action in the sector with Chern 
number C2 = 1. For SU{2) the corresponding action is equal to Sn'^/g'^, 
hence the tunnelling probability is enormously suppressed by the factor 
exp(— 47r/a,^) ~ 10~^™, where aw = fi'^/47r = a/sin^Ow — 1/29 {6w is 
the weak mixing angle). When the system is in the Higgs phase things are 
somewhat more complicated, since there is then no minimum for the Eu- 
clidean action with C2 = 1. Nevertheless, the tunnelling rate is suppressed 
by the same exponential factor. (For details, see [S2]-) 

Things are very different at high temperatures, corresponding to energy 
scales set by the height of the barrier. This is given by the static saddle 
point solution of the bosonic sector (Yang-Mills-Higgs equations), the so- 
called sphaleron solution. The energy of this spherically symmetric unstable 
solution is |3S! 

= W;(!!!£), (ni) 

aw rnw 

where the function varies from 1.56 to 2.72 as mH/mw varies from zero to 
large values. The height of the barrier is thus of order 10 TeV. At tempera- 
tures larger than the critical temperature Tc ~ mi^/^', 5- violating processes 
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may proceed through thermal activation over the barrier which is then much 
lower because the W boson mass eventually vanishes. The tunnelling sup- 
pression factor exp{—Esph/T) is no more effective in the symmetric phase. 
Detailed studies [SI] have established that the naive power counting estimate 
for the rate per unit volume in the unbroken phase 

T = 0{l){awT)^ (112) 

indeed holds. (For a review, see [22j.) We compare this with the expansion 
rate H of the Universe. During these early phases the energy density p is 
dominated by relativistic particles, thus p = (?*(7r^/30)T^, where (7* is the 
effective number of degrees of freedom per particle (2 for the photon). The 
Friedmann equation gives 

1 /2 

with gj = (9(10). For dimensional reasons, F should be compared with the 
time derivative of ( a typical particle density), i.e. with HT^. (We used 
that T is most of the time inversely proportional to the expansion factor.) 
The ratio 

r 4 Mpianck 

is much larger than 1 for Tc < T < 10^^ GeV. In the symmetric phase the 
5-violating reactions are thus fast enough to maintain thermal equilibrium. 

During this phase the B-\-L part of -B = ^{B — L) -\- -^{B -\- L) is therefore 
changed. At first sight one might expect that it is completely washed out, and 
that B = —L = \{B — L)initiai is established. Things are, however, somewhat 
more complicated, since only the fermion currents of the left-handed doublets 
have an anomaly. Moreover, electrical charge neutrality has to be taken into 
account. A more careful analysis shows |SS], that in good approximation 
one ends up above with 

B + L = --{B-L) (113) 

for the Standard Model with three generations and one Higgs doublet. (A 
mass correction jSH] changes this slightly.) Without a primordial B — L (e.g. 
from GUT physics) both B and L would vanish. These considerations are 
relevant for scenarios of leptogenesis (see [SE! ) ■ 
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Evolution below 

The electroweak phase transition cannot generate the cosmological baryon 
asymmetry, because the last of the following Sakharov conditions® is not 
fulfilled: 



• Violation of baryon number 

• Violation of C and CP 

• Deviations from thermodynamic equilibrium. 



The first two of these are obviously necessary. (The possibihty that the 

Universe is symmetric on scales much larger than the present horizon, but 
not in our visible part, will be briefly discussed in the next section) The third 
condition requires some explanation. A simple manner to understand its ne- 
cessity is this: In thermodynamic equilibrium one looses in a way the arrow 
of time (on cosmological short time scales), hence the guaranteed CTP sym- 
metry implies CP conservation. A more formal argument goes as follows: 
In equilibrium the density matrix that maximizes the entropy for given av- 
erage values of the absolutely conserved 'charges' Qi is the grand canonical 
ensemble 

Q^Z-^e-K^[-l3{H-Y.liiQi) 

i 

(Hi =chemical potentials for the conserved charges Qi). Suppose that B and 
L are independently violated, and that beside the electric charge Q no other 
charges like B — L are conserved. Because of electrical charge neutrality g 
the has to be the canonical ensemble: g — Z~^e~^^. But H is invariant 
under the CPT operation 9. Since B is odd under 9, we have 

{B) = Z-^Tr{Be-^") = Z-^Tr[9{Be-^")9-^] = -{B). 

So in thermodynamic equilibrium any previous baryonic and leptonic asym- 
metries would be washed out. Note that such an argument could not be 
used above, since not only B — L, but also all the anomaly free fermionic 
charges Nj, — ^B, i — e,ii,T (see (101) and (102)) are conserved. For this 
reason there remain, after imposing charge neutrality, three^ independent 

^In the words of Sakharov: "According to our hypothesis, the occurrence of C asym- 
metry is the consequence of violation oi CP invariance in the nonstationary expansion of 
the hot Universe during the superdense stage, as manifest in the difference between the 
partial probabilities of the charge conjugate reactions." 

^The differences Bi — Bj for the baryon numbers of the three generations are also 
anomaly free, however off-diagonal quark interactions violate all Bi (they preserve of 
course the total baryon number B). 
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chemical potentials. That B + L in (113) depends only on i? — L is the 
result of the approximations used in [33], j^H]- Similarly, baryonic and lep- 
tonic asymmetries may not be completely washed out if, beside the electric 
charge, there would exist other absolutely conserved charges whose chemical 
potentials are non-vanishing. Below Tc, as long as the 'sphaleron processes' 
are sufficiently fast to establish thermal equilibrium, the combination B + L 
is no more proportional to B — L. As a result, even ii B — L should vanish for 
all temperatures, a non- vanishing B would result (mass effects are important 
for this). 

In a strongly first order transition there would be a departure from ther- 
mal equilibrium and anomalous B and L violating processes could produce 
an asymmetry (with B — L = 0). Indeed, the transition would be very vi- 
olent through nucleation of bubbles in the new phase which would expand 
and collide. In order that the asymmetry is not erased after the phase tran- 
sition, the sphaleron energy at Tc must be sufficiently high. This leads to the 
following condition for the jump vt of the order parameter plotted in Fig. 7: 

vt/Tc > 1. (114) 

For a detailed discussion and references we refer to Sect. 7 of j^Z]. The results 
shown in Fig. 6 already imply that the electroweak transition in the Standard 
Model is too weak for baryogenesis. As discussed earlier, the transition even 
stops to be of first order for Higgs masses higher than ^ 72 GeV, much below 
the LEP limit mn > HO GeV. 

In summary, the baryon asymmetry of the Universe must have been pro- 
duced at energies much higher than the electroweak scale. Part of the B + L 
component will later be washed out by anomalous weak processes, but the 
B — L component would be preserved. Electroweak baryosynthesis is only 
possible in extensions of the Standard Model. The matter-antimatter asym- 
metry of the Universe is strong evidence for physics beyond the Standard 
Model. Possible scenarios have been reviewed many times; see, e.g., [TF] , 

m- 

3.3 QCD phase transition 

At temperatures of the order of the pion mass the pions are so abundant 
that they begin to overlap. As a result, the interactions of their constituents 
(quarks and gluons) become important. 

For a rough estimate, consider the number density of a free pion gas in 
the early Universe at T <^ rriN (nucleon mass). We can take their chemical 
potential equal to zero, because the pions are much more abundant than the 
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nucleons. Their number density is for T > m^: 



. _ 1 (27r)3 

The volume of a pion is (47r/3)r3, with the pion radius ~ 0.6 /m. We 
find rz^(T)t>7r = 1 at a temperature T ~ 285 MeV^, with an energy density of 
~ 850 MeV/fm\ 

On the basis of the asymptotic freedom of QCD one naturally expects that 
at high temperatures the hadronic matter behaves as an asymptotically free 
gas of quarks and gluons [quark-gluon plasma). Somewhere in the region 
T ~ there should be a phase transition or crossover between the low 
energy hadronic matter - with confinement and broken chiral symmetry - 
and the quark-gluon plasma. The study of the existence and nature of this 
phase transition is a central theme of QCD. One is, of course, also interested 
in the equation of state. 

Because of serious infrared singularities, perturbation theory is of limited 
value, even for the treatment of the asymptotic behavior of the quark-gluon 
plasma^. This is even more so when the phase transition is addressed. There- 
fore, the non-perturbative lattice approach is the only way to investigate the 
thermodynamics of hadrons from first principles. A lot of work has been 
invested in this domain over many years. Very useful general references are 
j7], [3H] and jHni- Since the critical temperature is close to the strange quark 
mass, it is important to simulate dynamical light quark flavors realistically. 
This has almost been achieved by now. A recent study of the MILC collabo- 
ration [40] comes to the conclusion that in the real world there is presumably 
no bona fide phase transition at the physical quark masses. In Fig. 8 we 
reproduce the result for the order parameter (ipip) for two degenerate light 
quarks (up and down) plus the strange quark with its physical mass. 

Fig. 9 shows the real part of the Wilson line or Polyakov loop, defined as 
(using the notation U^{x) for the link variables) 

Li^)=Tr\ li f/4(x,X4)l, (115) 
[2:4=1 J 

with Nt = {3 /a. This quantity is gauge invariant. The logarithm of its 
expectation value is proportional to the free energy of the system with a 
single heavy quark, measured relative to that in the absence of such a quark. 
(L) = can be interpreted as signalizing confinement. 

^For a brief discussion of the severe infrared problems, see Chapter 20 of Ref. Even 
resummed perturbation theory does not cure the infrared problems in QCD. 
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These and other results suggest that there is a crossover, rather than a 
sharp phase transition from the confined to the unconfined behavior. Similar 
results were obtained in jH]. 

Many cosmological speculations connected with the QCD transition (pos- 
sibility of inhomogeneous nucleosynthesis, strange quark nuggets, etc) have 
been discussed. Since it is not clear whether there are any observable foot- 
prints of this transition the Universe passed through (and for limitations of 
time), I leave this interesting subject. 

In the central regions of neutron stars matter is squeezed to densities 
several times higher than the nuclear density, pnuc = 2.8 x 10^"^ g cm~^ (see, 
e.g.. Chap. 6 of |42j). One expects that these regions are, at least for the most 
massive neutron stars, in the deconfined chirally symmetric phase. It is also 
expected that at relatively low temperatures the transition to this state of 
matter as a function of baryon chemical potential /ib (= '^l^'quark) is first order 
for realistic quark masses. It is difficult to perform lattice simulations for 
IXb 7^ 0, because the quark determinant (arising from the functional integral 
over the quark fields) becomes complex, and thus also the corresponding term 
in the effective action. For a clear discussion of these difficulties we refer to 
Sect. 5. 4. 3 of [Zj. Recently, significant progress has been made in algorithms, 
and realistic calculations may become feasible. 

3.4 Cosmic topological defects 

The formation of topological defects is well-known to condensed matter physi- 
cists. A classic example is the formation of vortices in type II superconduc- 
tors. As you all know these are described by the Ginzburg-Landau theory, 
which is a 'spontaneously broken' gauge theory (the Abelian Higgs model). 
The complex order parameter 0(x), which later turned out to be proportional 
to the gap function in BCS theory, plays in particle physics language the role 
of a Higgs field. I recall that the integer number in the fiux quantization is 
also the degree (winding number) of the map which associates to each direc- 
tion in the two-dimensional transverse space (topologically a circle 5"^) the 
asymptotic value of in the vacuum manifold. The latter is the circle on 
which |{/)| assumes the minimum of the Ginzburg-Landau (Higgs) potential. 

The analogous 't Hooft-Polyakov monopoles exist, for example as solu- 
tions of the SU (2)-Higgs model. The pole strength is then again topologically 
quantized, and the corresponding integer number is this time the winding 
number of the map which associates to each direction of 3-space (topolog- 
ically S"^) the asymptotic value of in the vacuum manifold, that is this 
time an 5*^. Note that the gauge group S'f/(2)acts on this vacuum manifold 
transitively and is thus topologically isomorphic to the homogeneous space 
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SU{2)/H, where H = U{1) is the stabihzer of any point of the vacuum man- 
ifold. For a general gauge theory with gauge group G, non-Abelian (regular) 
monopoles exist if the second homotopy group of G/H {H is again the sta- 
bilizer of the vacuum manifold) is non-trivial. Homotopy theory allows to 
compute these homotopy groups quite easily. 

Textures are also known in condensed matter physics. These are classified 
by a third homotopy group. 

The theory of topological defects in cosmology is extensively described in 
the monograph jlHI- On the basis of what we now know about the cosmic 
microwave background, the conclusion is inevitable that topological defects 
played no important role in large scale structure formation. This is the main 
reason why interest in the subject has declined. 



4 Vacuum energy problem and Dark Energy 

In the course of the various 'phase' transitions the Universe passed through, 
the free energy density changed by a huge amount in comparison to the 
present day cosmological energy density (which is close to the critical value). 
This gives me the opportunity to discuss the cosmological constant problem, 
a profound mystery indeed of present day physics. 



4.1 Vacuum energy and gravity 

When we consider the coupling to gravity, the vacuum energy density acts 
like a cosmological constant. In order to see this, first consider the vacuum 
expectation value of the energy-momentum tensor in Minkowski spacetime. 
Since the vacuum state is Lorentz invariant, this expectation value is an 
invariant symmetric tensor, hence proportional to the metric tensor. For a 
curved metric this is still the case, up to higher curvature terms: 

< Tfj^u >vac= dfiuPvac + higher curvature terms. (US) 

The effective cosmological constant, which controls the large scale behavior 
of the Universe, is given by 

A = S^Gp^ac + Ao, (117) 

where Aq is a bare cosmological constant in Einstein's field equations. 

We know from astronomical observations since a long time that p\ = 
A/SirG can not be larger than about the critical density: 

_ M 

Pcr^t - 
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1.88 X lO'^^hl g cm~^ 
(3 X 10-^eV)\ 



(118) 



where ho is the reduced Hubble parameter 



ho = ifo/(100 km s'^Mpc-^) 



(119) 



and is close to 0.7. 

It is a complete mystery as to why the two terms in (117) should almost 
exactly cancel. This is - more precisely stated - the famous A-problem. 

As far as I know, apart some unpublished remarks of Pauli in the early 
1920's jHj, the first who wondered about possible contributions of the vac- 
uum energy density to the cosmological constant was Zel'dovich . He dis- 
cussed this issue in two papers [45j during the third renaissance period of 
the A-term, but before the advent of spontaneously broken gauge theories. 
The following remark by him is particularly interesting. Even if one assumes 
completely ad hoc that the zero-point contributions to the vacuum energy 
density are exactly cancelled by a bare term, there still remain higher-order 
effects. In particular, gravitational interactions between the particles in the 
vacuum fluctuations are expected on dimensional grounds to lead to a gravi- 
tational self-energy density of order Gfj,^, where /i is some cut-off scale. Even 
for fi as low as 1 GeV (for no good reason) this is about 9 orders of magnitude 
larger than the observational bound. 

This illustrates that there is something profound that we do not under- 
stand at all, certainly not in quantum field theory ( so far also not in string 
theory). We are unable to calculate the vacuum energy density in quantum 
field theories, like the Standard Model of particle physics. But we can at- 
tempt to make what appear to be reasonable order-of-magnitude estimates 
for the various contributions. All expectations (some of which are dis- 
cussed below) are in gigantic conflict with the facts. Trying to arrange 
the cosmological constant to be zero is unnatural in a technical sense. It is 
like enforcing a particle to be massless, by fine-tuning the parameters of the 
theory when there is no symmetry principle which implies a vanishing mass. 
The vacuum energy density is unprotected from large quantum corrections. 
This problem is particularly severe in field theories with spontaneous sym- 
metry breaking. In such models there are usually several possible vacuum 
states with different energy densities. Furthermore, the energy density is 
determined by what is called the effective potential, and this is dynamically 
determined. Nobody can see any reason why the vacuum of the Standard 
Model we ended up as the Universe cooled, has - for particle physics stan- 
dards - an almost vanishing energy density. Most probably, we will only have 
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a satisfactory answer once we shall have a theory which successfully combines 
the concepts and laws of general relativity about gravity and spacetime struc- 
ture with those of quantum theory. 

4.2 Simple estimates of vacuum energy contributions 

If we take into account the contributions to the vacuum energy from vacuum 
fluctuations in the fields of the Standard Model up to the currently explored 
energy, i.e., about the electroweak scale Mp = Gp ~ 300 GeV {Op '■ 
Fermi coupling constant), we cannot expect an almost complete cancellation, 
because there is no symmetry principle in this energy range that could require 
this. The only symmetry principle which would imply this is supersymmetry, 
but supersymmetry is broken (if it is realized in nature). Hence we can at 
best expect a very imperfect cancellation below the electroweak scale, leaving 
a contribution of the order of Mp . (The contributions at higher energies 
may largely cancel if supersymmetry holds in the real world.) 

We would reasonably expect that the vacuum energy density is at least 
as large as the condensation energy density of the QCD transition to the 
confined phase of broken chiral symmetry. Already this is far too large, 
namely of the order ~ Aqj^-^/IOtt^ ~ 10~^GeV^, i.e., more than 40 orders 
of magnitude larger than Pcrit- Beside the formation of quark condensates 
< gg > in the QCD vacuum which break chirality, one also expects a gluon 
condensate < G^^Ga^u >~ ^qcd- This produces a significant vacuum energy 
density as a result of a dilatation anomaly: If denotes the "classical" trace 
of the energy-momentum tensor, we have 

TH = ©i: + ^GTGa,., (120) 

where the second term is the QCD piece of the trace anomaly {(^{g^) is the 
/3-function of QCD that determines the running of the strong coupling con- 
stant). I recall that this arises because a scale transformation is no more a 
symmetry if quantum corrections are included. Taking the vacuum expec- 
tation value of (120), we would again naively expect that < > is of the 
order Mp. Even if this should vanish for some unknown reason, the anoma- 
lous piece is cosmologically gigantic. The expectation value < G^^Ga^v > 
can be estimated with QCD sum rules [IZj, and gives 

< >««°™^ (350Me1/)^ (121) 

about 45 orders of magnitude larger than Pcrii- The expected energy density 
is larger than the energy-mass density in a neutron star. This reasoning 
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should show convincingly that the cosmological constant problem is indeed 
a profound one. (Note that there is some analogy with the (much milder) 
strong CP problem of QCD. However, in contrast to the A-problem, Peccei 
and Quinn ^HI have shown that in this case there is a way to resolve the 
conundrum. ) 

Let us also have a look at the Higgs condensate of the electroweak theory. 
Recall that in the Standard Model we have for the Higgs doublet $ in the 
broken phase for < $*$ >= the potential 

V^(0) = -\m'ct? + (122) 

Setting as usual (p = v + H, where v is the value of where V has its 
minimum, 



/2m2 



2-^/^G;p^/^ ~ 2466-6^, (123) 



V A 

we find that the Higgs mass is related to A by A = M'jj/v^. For (p = v we 
obtain the energy density of the Higgs condensate 

4 1 

V{<P = v) = -^ = -^MlMl = 0{M',). (124) 

We can, of course, add a constant Vq to the potential (122) such that it 
cancels the Higgs vacuum energy in the broken phase - including higher 
order corrections. This again requires an extreme fine tuning. A remainder 
of only 0(mg), say, would be catastrophic. This remark is also highly relevant 
for models of infiation and quintessence. 

In attempts beyond the Standard Model the vacuum energy problem so 
far remains, and often becomes even worse. For instance, in supergravity 
theories with spontaneously broken supersymmetry there is the following 
simple relation between the gravitino mass irtg and the vacuum energy density 

3 2 

Comparing this with eq.(30) we find 

Pvac 1 nl22 / ^9 \^ 



10122 (^y 



Pcrit ^TTlPl 

Even for rrig ~ 1 eV this ratio becomes 10^^. (m^ is related to the pa- 
rameter F characterizing the strength of the supersymmetry breaking by 
rUg = {AnG/Sy/'^F, so ~ 1 eV corresponds to F^/^ ^ lOO TeV.) 
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Also string theory has not yet offered convincing clues why the cosmo- 
logical constant is so extremely small. The main reason is that a low energy 
mechanism is required, and since supersymmetry is broken, one again expects 
a magnitude of order Mp, which is at least 50 orders of magnitude too large 
(see also jlHI)- However, non-supersymmetric physics in string theory is at 
the very beginning and workers in the field hope that further progress might 
eventually lead to an understanding of the cosmological constant problem. 

I hope I have convinced you, that there is something profound that we 
do not understand at all, certainly not in quantum field theory, but so far 
also not in string theory. 

For more on this, see e.g. 03] (and references therein). 

4.3 Astronomical evidence for Dark Energy 

A wide range of astronomical data support the following 'concordance' A- 
Cold- Dark- Matter (ACDM) model: The Universe is spatially fiat and dom- 
inated by vacuum energy density or an effective equivalent, and weakly in- 
teracting cold dark matter. Furthermore, the primordial fluctuations are 
adiabatic and nearly scale invariant, as predicted in simple inflationary mod- 
els. It is very likely that the present concordance model will survive. The 
evidence for a dark energy component of about 0.7 times the critical energy 
density is steadily increasing, and by the time the proceedings of this school 
will be available one can presumably make stronger statements. Readable 
reviews (for people outside the field) of the current evidence for a nearly ho- 
mogeneous energy density with negative pressure that dominates the energy 
content of the recent and future Universe are [20], Beside, there is a 

rapidly growing semi-popular literature on this exciting subject, confronting 
us with really profound fundamental physics questions. 
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